Modulational instability of spinor condensates 
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We demonstrate, analytically and numerically, that the ferromagnetic phase of the spinor Bose- 
Einstein condenstate may experience modulational instability of the ground state leading to a frag- 
mentation of the spin domains. Together with other nonlinear effects in the atomic optics of ultra- 
cold gases (such as coherent photoassociation and four-wave mixing) this effect provides one more 
analogy between coherent matter waves and light waves in nonlinear optics. 
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Recent experimental studies of Bose-Einstein conden- 
sation (BEC) in an optical trap have opened a new 
direction in the research of ultra-cold atomic clouds asso- 
ciated with their spin degree of freedom, or spinor BEC. 
An optical trap does not force atoms to align along the 
orientation of the strong confining magnetic field, as hap- 
pens in the case of a magnetic trap, allowing the study 
of atoms confined in all hyperfine states. Several recent 
theoretical studies have predicted a variety of novel phe- 
nomena that may occur in the spinor BEC, such as the 
propagation of spin waves and the existence of topologi- 
cal states -skyrmions, vortex states without a core 
In the case of spin-1 bosons such as ^^iVa, ^^K and Rb, 
the dynamics of the spinor BEC is described by the three 
spin degrees of freedom {mp = 1,0,-1 of the F = 1 
atomic hyperfine state) which are coupled parametrically. 
Depending on the parameters, such as scattering lengths, 
in an optical trap the ground state of the spinor BEC can 
be either ferromagnetic or antiferromagnetic ("polar"). 

In this paper, we demonstrate that the parametric 
coupling between the spin degrees of freedom provides 
a physical mechanism for the modulational instability of 
the ferromagnetic ground state, for large enough densities 
of spinor BEC. In contrast, the antiferromagnetic state 
is always modulationally stable. This effect is reminis- 
cent of the quasiparticle instabilities in two-component 
homogeneous condensates that are known to occur only 
for certain ratios of the inter- and intra-component inter- 
action strengths It also suggests one more example 
of a deep analogy between the coherent matter waves 
and light waves in nonlinear optics, along with the al- 
ready studied cases of four-wave mixing [ |lO| and atomic- 
molecular photoassociation [ pd[ |. 

Model. We consider an atomic spinor BEC in an opti- 



cal trap in the magnetic- field-free case. The Hamiltonian 
for the spinor BEC in the optical trap has the form H: 
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' {v)'ip^b^{v). In the equations above 
7i^V^/2TO -I- Vt{v) is the single-atom Hamil- 



tonian, and Vr(r) is the trapping potential created by 
an optical field. Without loss of generality, we restrict 
ourselves to the case of hyperfine atomic state F — 1 
with the corresponding sub- levels: rap — —1, 0, 1. 

The nonlinear interaction between different spin com- 
ponents of the condensate is governed by the spin- 
independent and spin-dependent interaction strengths, 
Co = 47r?i^(2a2 + ao)/3?7i and C2 — 47r?i^(a2 — ao)/3m, 
respectively. Here ag and 02 are the s-wave scattering 
lengths for two colliding atoms with total angular mo- 
mentum Ftot = and -Ftot = 2. The coefficient C2 is 
calculated to be positive for the polar spinor condensate 
(e.g., "^^Na) and negative for a ferromagnetic condensate 
(e.g., ^'^Rb) (see, e.g., g). Whether the spinor conden- 
sate is ferromagnetic or antiferromagnetic arises from the 
fundamental gauge symmetry of the system, which in 
turn is dependent on the sign of the coefffcient C2 (cq 
is always positive for condensates with repulsive interac- 
tions, such as the "^^Na and considered here). 

From Eq. (|l]) we can derive three coupled field equa- 
tions for the case under consideration: 



"^^-gf = ^'^1 + C2(^l'0l + V'dV'o - '0^^i'0-l)'0l + C2'0ll'0O: 

^ ^^^^ C2(^/;liV_i + V^J^o - + C2'lp\'lpl, 
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where L = [-^V^ + VT{r)\ + co(V'LiV'-i + V-JV-o + 

Next, we introduce the foUowing Hnear field 
transformation for the quantum-field components, and 
treat the components as the corresponding mean fields: 
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As will become clear below, the advantage of this basis 
is that the stationary solutions for all three spinor com- 
ponents can be described by a single equation. 

The coupled dynamical equations in the new basis have 
the following form: 
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where this time C = [-jV^ + Vrir)] -I- cq ( 1 0+ P -I- 1 00 P + 
|0_P). The wave functions, time, spatial coordinates, 
and interaction strengths are measured in the units of 
(?i/TOtj)-3/4, uj'\ (?i/mw)i/2, and ihuj)-^{h/muj)-^^^, 
respectively, where lo is the axial trapping frequency. 

Stationary stoics. The dynamics of the spinor conden- 
sate described by Eqs. (Q) is, in general, spin mixing. 
However, any stable stationary solution of the system 
represents a non-spin-mixing, or spin-polarized state 
of the system. Such stationary states have a constant 
population of each spin component and can be found by 
introducing the following ansatz: 
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where j — (-f,— ,0), 9± are the relative phases of each 
component with respect to 4>o{0q = 0), and iij are the re- 
spective chemical potentials. Finally, + rig + n_ = n 
is the total density of the spinor condensate. 

Upon substitution of the ansatz (||) into Eqs. (||), it 
becomes apparent that the stationary solutions can ex- 
ist in this dynamical system only under the condition 
yLt_l_ = ^_ = /j,o = Ai- Another condition is that the spin 
components of the condensate are locked in phase. For a 
stationary state to exist, the variables 29- and 26'+ can 
take only two distinct values, or tt. In such a station- 
ary state, the spinor eigenfunctions y^rij are eigenmodes 
of the same effective potential created both by the optical 
trap and by the nonlinear interaction of all spin compo- 
nents. Moreover, all three eigenmodes correspond to the 
same eigenvalue, /i. This means that these eigenfunc- 
tions are proportional to each other and therefore can be 



e'"+ = 1; r+ = ro + = 1/2, and 
= 1/2, and 
= 1/2, and 



presented in the form: nj{r) = rjn{r), where constant co- 
efficients rj represent the population of each spinor com- 
ponent in a steady state, with ?'+ + r_ + rp = 1. The 
spatial profiles of all three stationary spinor components 
obey the same time-independent equation: 

(^~^V^-^i + VT{r)^V^ + Xn{r)V^ = 0. (6) 

Depending on the relative phases, there exist four dif- 
ferent phase-locked steady state solutions of Eqs. @) with 
different populations in each spin component and differ- 
ent values of the coefficient x in Eq. (g|): 

Case 1. e'^" — 1, e*^+ — i. tq and r± are constrained 
by r+ + r_ -I- tq = 1, and x — cq- 

Case 2. e^^- =1, 

X = Co+C2. 

Case 3. e'^- = i, 6*"+ = i; r_ = + r^ 

X = Co C2. 

Case 4. e'^- = i, e'^+ = 1; tq = r_ + 

X = Co -t- C2. 

Modulational stability analysis. The stationary so- 
lutions described above may correspond to different 
metastable states of the spinor condensate, provided that 
they are linearly stable. Linear instability, i.e. exponen- 
tial growth of the modulation of the stationary homo- 
geneous condensate, was previously hinted to be respon- 
sible for complex spatial modulations of the condensate 
in a trap, that ultimately lead to the destruction of the 
non-spin-mixing state However, no stability analy- 
sis of the spinor EEC has been carried out previously. 
Similar phenomenon occurring due to the nonlinear in- 
teraction of light beams (and pulses) in optical media is 
called modulational instability (MI) and is well studied in 
the context of nonlinear optical fibers ||l^ . 

To perform the MI analysis for the spinor BEC, we 
first note that the stationary homogeneous (or constant 
density) solutions corresponding to the phase relations 
in the cases 1-4 above have the form: = n/x- Next, 
we add a small (generally complex) perturbation to the 
homogeneous solutions, taking the functions 4>q and (f)± 
in the form: 



0j(r,<) = (./n^ + (50,)e 
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Substituting Eqs. (0) into Eqs. (jj), omitting the 
terms containing Vr, and linearizing around the homo- 
geneous solutions, we obtain the dynamical equations 
for the perturbations S(j)j (r, t) , which can be used to 
analyze the stability of the homogeneous solutions in 
the cases 1-4 against the growth of periodic pertur- 
bations. If the perturbations are taken in the form 
S(j)j — (uj + jwj)cos(kr)e"*, where k = {kx,ky,kz), 
these equations become: Afl^ — 0, where fl — 
(w+, w+, Uq, Vq, u_, «_), the matrix A being too 
cumbersome to write out here. If the characteristic equa- 
tion, det A = 0, resolved with respect to the perturbation 
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frequency w, has a real or complex root for some real pos- 
itive k^, the spinor condensate is modulationally unsta- 
ble. In general, such an equation is of the sixth order in 
(jO. For simplicity, we can assume one of the populations 
rj equal to zero, and the other two equal to each other. 
Then for case 1, all possible eigenvalues are given by: 

2 V 2 Coy 

From these expressions, and keeping in mind that 
|c2/co| < 1, one can see that the real positive values 
of up , and hence the MI of this solution, can occur only 
for the ferromagnetic state, i.e. when C2/C0 < 0. 
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FIG. 1. Time evolution of the central density n(0) and 
atom number, N, of the components of the ID spinor conden- 
sate (p- (dots) (fi+ (dashed), (po (dot-dashed) and the total 
density and atom number (solid). Parameters are e'^~ — 1, 
e'*+ =i,ro = 1/2, r± = 1/4, xo = 0.0528, X2 = -0.00048, 
and a = 1; <j!>_|t=o ~ <^+|t=o- The total number of atoms in 
the system is A'^ = 5 x 10^ corresponding to the ID chemical 
potential fi ~ 125. 

Carrying out an identical analysis for the cases 2-4, we 
find the following possible eigenvalues: 



0.,^ = -2(^^-2^^)', (9) 

V 2 Co + C2 / 

u;i = -2 (^y + 2^^ . 

Here all w| are negative and therefore, MI does not occur 
in any of the cases 2-4, neither for a ferromagnetic nor 
for a polar state. 

It is possible to show that the conclusions of the MI 
analysis above hold in the most general case, when none 



of the populations rj are zero. Thus homogeneous solu- 
tions in the polar state never experience MI. 
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FIG. 2. Time evolution of the ID (/f>-(z) component for 
the case 1, showing the spatial development of the MI. The 
(l)+fi(z) components behave similarly. Parameters are as in 
Fig. 0. 

Numerical simulations. The MI analysis for the ho- 
mogeneous condensate does not, strictly speaking, apply 
to the trapped condensate. However, it can serve as an 
indication of the spinor condensate behavior since the in- 
stability of the homogeneous condensate is bound to trig- 
ger the formation of the complex patterns in the trapped 
condensate cloud if the characteristic spatial extend of 
the condensate, I (in dimensionless units), is larger than 
the largest length scale of the spatial patterns due to 
the MI, L — kz}„. For the only case when the MI of 
the spinor system does occur, this condition becomes 

I > (m|c2/co|)-1/2. 

To test the results of our stability analysis and to 
demonstrate the effect of the MI on the trapped spinor 
condensate, we carry out numerical simulations of the dy- 
namical equations (^ . Because of the anisotropic nature 
of the optical trapping potential, the cigar shaped EEC 
is assumed to be quasi-one-dimensional and hence we use 
the ansatz (pj{r, t) = ^{x, y)(j)j{z, t), where z is the direc- 
tion of weak confinement, and '^{x,y) is the the wave- 
function of the two-dimensional harmonic oscillator. This 
leads us to the one-dimensional (ID) dynamical spinor 
system for (t)j{z,t) which is identical to the system (0) 
with £ = i {-d^/dz^ + z^) + xo(|0+P + \<l>o? + 
The dimensionless interaction coefficients are xo,2 = 
Co, 2a, and a = ^ \^{x,y)\'^dxdy / ^ \^ {x , y)]"^ dxdy is the 
transverse structure factor. Using the initial conditions 
specified by: 

ct^A^) = ^jr,n{zV''^ (10) 

where n[z) is the ID spatial profile determined from Eq. 
(^), we solve the ID equivalent of Eqs. (^) numerically 
for the cases 1-4 of the stationary phase- locked solutions. 
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In agreement with the analytical results, we only observe 
MI for the ferromagnetic state [02/ cq < 0) in the case 
1. The results of a representative calculation for the fer- 
romagnetic state (corresponding to ^^i?6) are shown in 
Figs. 1^ and |^. In all other cases the spinor BEC remains 
stable to periodic modulation of its components. 

Our results show that the effect of the MI on the 
trapped spinor BEC is twofold. First, the instability 
destroys the non-spin-mixing stationary state, leading 
to population transfer between the spin components, as 
shown in Fig. (|^). Second, the MI causes periodic spatial 
modulation of the condensate which, being confined by 
the trap, grows to be chaotic with time (see Fig. ^). In 
the short term, MI causes the spatial fragmentation of 
the spin domains shown in Fig. |3| for the original spinor 
components, ipj. 
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FIG. 3. Spatial intensity profiles of the original spinor 
components, ipj, demonstrating spin domain fragmentation. 
V'-(z) (dots), tp+{z) (dashed), Vo(z) (dot-dashed) and the 
total density (solid). Parameters are as in Fig. |l| 

Finally, to check that our results are dimensionally in- 
dependent we have also numerically analyzed the spinor 
system in two dimensions, and again found that the MI 
only occurs in the ferromagnetic state for the stationary 
solution of case 1, as shown in Fig. ^. 

In conclusion, we have predicted analytically and 
demonstrated numerically the possibility of MI of the fer- 
romagnetic ground state of the spinor BEC. This effect 
resembles the parametric MI in birefringent optical fibers 
with the Kerr nonlinearity, and it provides one more di- 
rect analogy between matter wave, and nonlinear, optics. 
Our results show that the ferromagnetic spinor BEC is 
an ideal candidate for the first experimental observation 
of MI, spin domain fragmentation, and spatio-temporal 
chaotic dynamics of matter waves. 

The authors gratefully acknowledge a reading of this 
manuscript by Dr. H. Pu. This work has been partially 
supported by the Performance and Planning Fund of the 
Institute of Advanced Studies. 
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FIG. 4. Developement of MI in the 2D case, shown for the 
spatial intensity profile of the spin component. Intial con- 
ditions are calculated from Eq. (^. Parameters are as Fig. ^ 
with the structure factor a^n ~ 1. 
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